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Abstract. We study the non-perturbative phenomena of Dynamical Mass Generation and 
Confinement by truncating at the non-perturbative level the Schwinger-Dyson equations in 
Maxwell-Chern-Simons planar quantum electrodynamics. We obtain numerical solutions for the 
fermion propagator in Landau gauge within the so-called rainbow approximation. A comparison 
with the ordinary theory without the Chern-Simons term is presented. 



1. Motivation 

Dynamical Mass Generation (DMG) and Confinement are two emergent plienomena of Quantum 
Chromodynamics (QCD) responsible for the nature of the hadronic spectrum. These can be 
studied by means of the corresponding Schwinger-Dyson Equations (SDE), as an alternative to 
the formulation of QCD on the lattice. The SDEs form an infinite tower of relations among the 
Green's functions of a given quantum field theory, whose nature is genuinely non-perturbative. 
We choose to study the SDE in a simplified model that shares with QCD the phenomena of 
Confinement and DMG, namely, quantum electrodynamics in the plane, or QED3. Applications 
of QED3 also extend to condensed matter systems, including high-Tc superconductivity, the 
quantum Hall effect and, very recently, graphene. Moreover, it is interesting to include the 
Chern-Simons (CS) term for the gauge field, which explicitly breaks parity and induces a 
topological mass for the photons. 

Incorporating the CS term in studies of DMG in QED3 have found that when the 
CS coefficient reaches a critical value, chiral symmetry is restored in a first order phase 
transition [Il[2]. In this work, we are interested in the effect of the CS term regarding confinement 
in quenched QED3. In Sec. 2 we present our model with a reducible representation of (2-1-1)- 
dimensional quantum electrodynamics with a CS term. Solutions to the gap equation in the 
rainbow approximation are presented in Sec. 3. The question of confinement in QED3 with a 
CS term is addressed in Sec. 4. Finally, in Sec. 5, we present our conclusions. 



2. Model 

We study quantum electrodynamics in (2+l)-dimensions. For fermions, there exist two 
irreducible representations of the Euclidean Clifford algebra {7^,71/} = 25^^. In any of 
these representations, it is impossible to define chiral symmetry [31 [H [5]. Furthermore, the 
fermion mass term, regardless its origin, is parity V non invariant. One still can construct 
a parity-preserving Lagrangian considering two different species with a relative sign between 
their masses [6l [7|. As a result, two chiral transformations can be defined [7]. These two 
species are conveniently merged into a single four-component spinor, making use of a reducible 
representation of the Dirac matrices. 

As compared with its four-dimensional counterpart, only three Dirac matrices are required 
to describe the dynamics of planar fermions. Once we have selected a set of matrices to write 
down the Dirac equation, say {70,71,72}, two anti-commuting gamma matrices, namely, 73 and 
75 remain unused. Hence the corresponding massless Dirac Lagrangian is invariant under the 
chiral-like transformations ~^ e^'^'^^il) and ■0 — ^ e^^'^^ip, that is, it is invariant under a global 
U{2) symmetry with generators 1, 73, 75 and [73, 75], corresponding to the interchange of fermion 
species. This symmetry is broken by an ordinary mass term rrieiptp. The order parameter of the 
symmetry breaking is the condensate {ipip)e = {0\ipip\0). There exists a second mass term which 
is invariant under the "chiral" transformations, referred to as the Haldane mass term nioipTip [8] 
with r = [73,75]/2 = diag{I,—I), which is associated with the condensate {4''ip)o = (O|V'T'0|O). 
The ordinary mass term is even under parity V transformations, but the Haldane mass term is 
not. We use the subscripts e for even and o for odd in quantities of the model. The Haldane 
mass term radiatively induces a parity-odd contribution into the vacuum polarization, which 
can be traced back to an induced Chern-Simons interaction 

iff 

C-CS = -^£:fiup\F^p ■ (1) 

Such a term is parity non invariant, and despite the fact that it is not manifestly gauge invariant, 
the corresponding action respects gauge symmetry. The parameter 9 induces a topological mass 
for the photons. Thus, the model we consider is described by the Lagrangian 

C = ^{i$ + eA + me + Tmo)^ + \Ff,yF^y + ^{d^A^f -^-^e^^ypApFup. (2) 

4 Zi^ 4 

There are many condensed matter systems which can be described by this Lagrangian, for 
which the physical origin of the masses depends on the underlying system, including high-Tc 
superconductivity, the quantum Hall effect and, very recently, graphene. 

Written in this form, neither m-g nor ttt-o correspond to poles of the fermion propagator. We 
introduce the chiral projectors x± = (1 =t t)/2 which have the properties x± = X± iX+X- = 
, x+ + X- = 1- These allow us to define the right-handed ■0+ and left-handed ip- fermion fields 
as '0-1- = XitV'- The x± project the upper and lower two component spinors (fermion species) 
out of the four-component "0, such that the fermion sector of the Lagrangian ^ can be cast in 
the form 

Cf = i^+i^ - 'm+)i>+ + - ■m_)'il^_, (3) 

with m± = TTig iruo- This Lagrangian describes two fermion species of physical masses and 
m_, respectively. These masses break chiral symmetry and parity at the same time. Moreover, 
the effect of the parity-violating mass is seen to remove the mass degeneracy between the two 
species. In what follows, we will be interested in the analytical properties of the dynamically 
generated fermion propagator associated with these physical masses. The analytical structure 
of the propagator can be studied from the corresponding Schwinger-Dyson equation, which 



corresponds to the expression 



where T^{k,p) and A^^{k — p) are, respectively, the full fermion-photon vertex and the full 

photon propagator, which obey themselves their own SDE, and Sp{p) and Sp\p) stand for the 
full and bare fermion propagator. In QED3, the coupling has mass-dimension one. Moreover, 
since the theory is super-renormalizable, becomes the natural scale of massless QED3, which 
is directly connected to the scales of confinement and dynamical chiral symmetry breaking. In 
this work, we write the relevant mass scales in units of = 1. 

From the Lagrangian ([2]), we observe that inverse fermion propagator has the form 

Sp\p) = Ae{p)f + Ao{p)Tf - Be{p) - Bo{p)t . (5) 

The scalar functions Ae^o{p) and Bg^oip) can be expressed in terms of the fermion wavefunction 
renormalizations Fe.o{p) and the mass function Mf>^o{p) as Ae^oip) = l/-^e,o(p) and i3e,o(p) = 
Me,oip) / Fe^oip) both in the even and odd sectors. The bare propagator corresponds to the 
values A^\p) = 1, a'^\p) = 0, B^\p) = me, B^\p) = ruo- Rather than working with 
parity eigenstates, we find it convenient to work with the chiral Lagrangian The chiral 
decomposition of the fermion propagator becomes 

A+{p)p + B+{p) A^{p)f + B_{p) 

Al{p)p^ + Bl{p)^+ Al{p)p^ + B^_{p)^- 



+ al{p)]x+ - K (pV + ^ , (6) 



A(o)fol - i (5 - ^] - ^""P'^P^ (7) 



where K± = it Kq for all relevant even and odd quantities of the model. The inverse 
transformations are simply = {K^ + K^)/2 and Kq = — K-)/2. In this basis, the 

bare quantities are Af{p) = 1 and Bfip) = m±. On the other hand, the photon propagator 
associated with the Lagrangian ([2]) in the Landau gauge, ^ = 0, takes the form 

+ 02 l^V ^2 ; ^2(^2 + 
3. Gap Equation 

We consider the rainbow-ladder truncation scheme, in which one makes the replacements 
and A^^ ^ A^). 

With the chiral decomposition of the fermion propagator, Eq. ([6]), both in the left- and right- 
handed sector, the SDE corresponds to a matrix equation which can be converted into a system 
of scalar equations for A±{p) and B±{p) after multiplying it by 1 and p and taking the trace. 
In the limit m± = 0, this leads to the decoupled system of equations 

A^{p) = l + ^fd^kalik) tf''j} T4-2 fd'kaiik)- ' 



with Q = l/(47r) as usual. Because of the effects of the CS term, the A±{p) functions become 
not trivial. Thus, in covariant gauges these never decouple from the equations for B±[p) as long 
as 7^ 0. We numerically solve the above systems of equations in each sector varying 9. 
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Figure 1. Functions A±{p) for various values of 9. 
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Figure 2. Functions B±{p) for various values of 9. 



In Fig. [T] we show the functions A±{p) for different values of 9. The small deviation from 
the tree-level value for small momentum in each case is due to the CS term. As 9 increases, 
both functions exhibit a plateau for small p and smoothly tend to 1 in the ultraviolet. The 
deviation is more pronounced in A^{p). Figure [2] shows the result for B±{p). We observe that 
the height of the plateau of B^{p) at small momentum decreases from the parity-even result as 
9 gets bigger, but the height of the plateau of B^{p) exhibits the opposite behavior. The effect 
is more visible in B^{p). In fact, there is a critical value 9^ — 8 x 10^^ above which B^{p) 
strongly changes its behavior and the plateau becomes negative, as shown in Fig. [3l 

This sudden drop has strong implications regarding chiral symmetry restoration in the 
model [HE]- In order to illustrate this, in Fig. [U we show the mass functions M-i-{p) = 
B±{p)/A±{p) for different values of 9. At 9c, the plateau of M+(p) also drops to negative 
values, whereas for M_(p), it continues to increase with 9. To ensure that these results are 
physically sensitive, we study the asymptotic behavior of M-i-{p). In the infrared, the height of 
the plateau can be considered as an order parameter for dynamical chiral symmetry breaking. In 
this connection, in Fig. [5] we draw the dependence of = M-i-(O) as a function of 9 below and 
above criticality. This parameter can be regarded as the dynamical mass of the corresponding 
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Figure 3. Critical behavior of B^{p) for ~ 8 x 10 
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Figure 4. Mass functions M±{p) for various values of 9. 



fermion species. We observe that the role of the CS coefficient is to remove the mass degeneracy 
between fermion species as long as 9 < 9c- There is a light and a heavy species. At 9c, 
however, there is a drastic change in this behavior, the light species develops a negative mass, 
which in absolute value is the same as its heavy cousin, as can be appreciated in the dotted 
curve in the same graph. This implies that above this 9c, the would-be parity preserving mass 
He = + /i-)/2 vanishes, and hence chiral symmetry is restored [HH]. Nevertheless, we want 
to emphasize that in this model, //g does not correspond to a pole in a propagator and hence to 
a physical mass. Physical masses iJ,± are generated for arbitrarily large values of 9. 

The discontinuous behavior of the dynamically generated masses suggest a first order phase 
transition for the dynamical breaking of chiral symmetry. Whether there is chiral symmetry 
restoration can better be seen from the behavior of the condensates, which are the true order 
parameters. These can be extracted from the solution to the SDEs as 

< >±= / dkk^(ji{k) . (9) 

Jo 

The results are shown in Fig. [6l Notice that in the parity preserving case, 9 = 0, the value of 
the condensate for each species is half the value of the ordinary condensate [9l[T0]. The reason is 
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Figure 5. /ij- as a function of 6. 
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Figure 6. Chiral condensates as a 
function of 9. 
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Figure 7. Comparison of the dynamically generated fermion masses for the full system ([8]) and 
the perturbatively reduced system p^ . 



that when fermion species are degenerate in mass, each one of them, corresponding to a fermion 
pairing mode ipiip^ or ipfipi (arrows indicate spin orientation), contributes one half to the parity 
preserving condensate. 

Finally, let us consider wavefunction renormalization effects. Regarding the wavefunction 
renormalization in the dynamically generated mass functions M±(p), we compare the solutions 
of the full systems in the previous subsection, Eq. ([8]), against the solution for the simplified 
equations resulting from the perturbation theory motivated ansatz = 1, namely 

a /■ 3 B±{k) 1 f jSi, ^ ' q) 

The behavior of the solutions in this truncation is qualitatively similar to that of the complete 
system ([8]). The critical behavior of the dynamical masses is not washed out by the effect of 
wavefunction renormalization, as can be seen in Fig. [71 where a comparison of fi± and 
between the full result and the perturbatively motivated one is shown. The behavior is virtually 
the same in both truncation schemes. The critical value 0p^'^^ is of the same order of 9c and thus 
the effects of A±{p) ~ 1 are not significant. This means that even with the inclusion of the CS 
term, the rainbow-ladder truncation is still reliable in Landau gauge. 
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Figure 8. Confinement tests for the spatially averaged Schwinger functions. 



4. Confinement 

In ordinary quenched QED3, there is confinement. Here we are interested in the extent up 
to which this scenario persist in the presence of a CS term. Whether a solution to the SDE 
supports confinement can be tested by means of the violation of the Osterwalder-Schrader axiom 
of reflection positivity, which states that the spatially averaged Schwinger functions 



A±(t) 



(fx 



(2^ 



'a. 



(11) 



should be positive definite if they are related to a stable asymptotic state. In our case, we 
construct the functions A±(t) inserting the solutions for the SDE in the above expression, upon 
which we perform a confinement test. In Fig. [8] we plot the logarithm of the absolute value 
of each of these functions. The oscillatory behavior inferred from the pronounced peaks in the 
graphs reveals that even in the presence of a CS term, quenched QED3 exhibits confinement. 
This behavior corresponds to two complex conjugate mass like singularities with complex masses 
m = azizib for each species, which fit to 



A±(t) ~ C±e-''^^cos{b±t + 6±) . 



(12) 



Thus, confinement is observed in our model. 

Following Ref.flT], we perform a similar analysis on the vector parts of the propagator, a^. 
We define new set of spatially averaged Schwinger functions, 



(2^ 



„i(tpo+2:-p) V 



(13) 



and insert again the solutions found in previous sections. In Fig. [9] we draw the logarithm of 
the absolute values of Q.±{t). Again, the oscillations hint confinement. Moreover, we obsverve 
that ^±{t) can be fitted according to 



(14) 



which also unveils the complex pole mass structure of the propagator. 
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Figure 9. Spatially averaged Schwinger functions involving the vector parts of the fermion 
propagator. 



5. Conclusions 

We have studied the dynamical generation of masses and confinement in Maxwell-Chern-Simons 
QED3. There exists two types of fermions within the four-component spinor formalism. These 
species are non-degenerated in mass. The origin of these physical masses is two-fold: On the 
one hand there is a parity-preserving contribution mg coming from dynamical chiral symmetry 
breaking, on the other hand, there is the CS-induced Haldane mass term ruo- We have explored 
the consequences of confinement in our model. 

For 9 > 6c the masses for each species are equal in magnitude, fi+ = — and we do not 
generate the ordinary mass fie = (/^+ + /^~)/2 that preserves parity. This is the reason why 
there is restoration of chiral symmetry and we only generated masses induced by the CS term. 
The phase transition for the mass is of first order. Regarding confinement, according to 
tests performed with and o"^, we observe that quenched QED3 with a CS term also shows 
this feature, regardless the value of the CS coefficient. This is in contrast to the case in which 
vacuum polarization effects are considered in the leading order of the approximation 
where there is deconfinement for arbitrary finite values of 9. 
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